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Aims. The identities of the main processes triggering and quenching star-formation in galaxies remain unclear. A key stage in evolution, however, 
appears to be represented by post-starburst galaxies. To investigate the prevalence of these galaxies and their impact on galaxy evolution in general, 
S^/ ■ we initiated a multiwavelength study of galaxies with k + a spectral features in the well-studied COSMOS field. 

Methods. We examine a mass-selected sample of k + a galaxies in the COSMOS field at z = 0.48- 1.2 using the spectroscopic zCOSMOS sample. 
To classify galaxies in their k + a phase, we use a spectroscopic criterion, based on the amplitude of the US absorption line and the absence of 
the [Oil] emission line. We develop our analysis for a well-defined sample by imposing stringent confidence levels on the spectroscopic redshifts 
Q ■ and spectral measurements. We compare our results for two mass-selected samples of star-forming and quiescent galaxies selected using a purely 
S—j ' spectral classification scheme from the 10,000 zCOSMOS catalogue (i.e. , based on measurements of 4000 A break and EW[OII]). 
CZ) , Results. In our mass-limited sample, k + a galaxies occupy the brightest tail of the luminosity distribution. They are as massive as quiescent 
, galaxies and populate the green valley in the colour versus luminosity (or stellar mass) distribution. A small percentage (< 8%) of these galaxies 
have radio and/or X-ray counterparts (implying an upper limit to the SFR of ~ 8A4 yr~' ). Over the entire redshift range explored, the class of post- 
starburst galaxies is morphologically a heterogeneous population with a similar incidence of bulge-dominated and disky galaxies. This distribution 
does not vary with the strength of the Hi5 absorption line but instead with stellar mass in a way reminiscent of the well-known mass-morphology 
relation. The results about the incidence of asymmetries and the concentration of the light distribution derived from HST/ACS images imply that 
this galaxy population possibly represents an intermediate stage of galaxy evolution. Although k + a galaxies are also found in underdense regions, 
they appear to reside typically in a similarly rich environment as quiescent galaxies on a physical scale of ~ 2-8 Mpc, and in groups they 
show a morphological early-to-late type ratio similar to the quiescent galaxy class. With the current data set, we do not find evidence of statistical 
significant evolution in either the number/mass density of k + a galaxies at intermediate redshift with respect to the local values, or the spectral 
Q\ properties, although more solid results on this and other aspects will be obtained following the completion of the survey. 

f^*) ' Conclusions. Several mechanisms related and unrelated to the environment are at work in quenching star-formation activity in galaxies on short 
Q\ | timescales (< 1 Gyr). Those galaxies, which are affected by a sudden quenching of their star-formation activity, may increase the stellar mass of 
f^**) , the red-sequence by up to a non-negligible level of ~ 10%. 

^> , Key words. Galaxies: formation - galaxies: evolution - galaxies: fundamental parameters - galaxies: mass function - cosmology: observations 

X: 

C$ ' 1. Introduction iBrinchmann & Ellisl I20 00). observational evidence exists that 

many galaxies have e xperienced morph ological transformations 

Numerous studies have demonstrated that a bimodal popula- at various re dshifts. Brinchmann & Ellis (2000) reported that 

tion of galaxies existsupto high redshift (e.g., | Franzetti et al. | a decline has occurre d in the mass density of irregular galax- 

|2007t | Mignolietal. | |2009|). The first population consists of iefj between z ^ f and to day, in contrast t0 an increase for 

galaxies that are actively star-forming, gas-rich, rotationally- regular galaxies . Further support of this redistribution between 

supported, and morphologically disk-dominated. The second galaxy types has been provided by study of th e stellar mass 



population represents typically quiescent galaxies that are function by many di fferent surveys, e .g., K20 (Fontana et al. 
gas-poor, pressure-supported, and have spheroidal morpholo- COM BO-17 (iBorch et al] ITO)" D EEP2 ( Bundy_etaL 

gies. As noted by many studies (e.g., [Bower et al. | LL998fc |2Q06h. SWIR E- VVDS-CFHTLS (lArnouts et all 120071) . SPSS 



(Martin et al. 2007), and VVDS (Vergani et al. 2008). It is im- 



Send offprint requests to: D. Vergani, e-mail: por tant to assess the role of environment in these transforma- 

daniela . verganiOoabo . inaf . it tions of galaxy type For i nstance , it has been shown that galaxies 

* Based on data obtained with the European Southern Observatory ^ dense [Qm haye modifl ed distributions of HI, the primary 
Very Large Telescope, Paranal, Chile, program 175.A-0839. 
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source of fuel for star-formation, and truncated star-formation 
activity (e.g.. iGiovanelli & Havnesl fl98l ICavatte et al.1 Il990t 
Koopmann & Kennev 2004). However, several other mecha- 
nisms may also contribute (see below). Thus, important ques- 
tions in galaxy formation and evolution are whether star-forming 
galaxies are the precursors of the quiescent galaxy population 
observed at the present day, or whether there are alternative paths 
in galaxy formation, and, if so, what are the primary processes 
involved and how are they controlled? 

The main limiting factor in answering this question is that we 
still have no true understanding of what triggers star formation 
in galaxies and why this activity is suppressed in some of them. 
There is a class of galaxies that might shed light on this issue. 
They are called k+a or E+A galaxies (depending on the author), 
or simply galaxies in the post-starburst phase. In this work, we 
adopt the spectral terminology k+a instead of the morphological 
E+A nomenclature. 

Post-starburst galaxies, discovered by Dressier and Gunn 
(1983), exhibit peculiar spectral characteristics, i.e., strong 
Balmer absorption lines, indicative of an intense star-formation 
epoch in the past billion years, and an absence of emission lines, 
a signature of a lack of ongoing star-formation. On the basis of 
these spectral signatures, they have been assumed to represent 
galaxies in a transitional stage between being a star-forming, 
late-type galaxy and a passive, early-type system. 

Various mechanisms have been proposed to explain the char- 
acteristic spectral features of post-starburst galaxies. Cluster- 
related mechanisms, su ch as ram-pressure gas stripping, harass- 
ment, or strangulation dGunn & Gottlll972t lLarson etal 1980; 
Balogh & Morris 2000) are efficient in suppressing the star- 
formation activit y of a galaxy. Events such as galaxy mergers 
and i nteractions (Toomre & Too mrelll972l; iBarnes & Hern quist 
|1992|) . commonly found in both the field and in clusters (al- 
though more efficient in the field), have been proposed. The 
large fraction of morphological irregularities (of similar phys- 
ical timescale as a post-starburst phase) and galaxy compan- 
ions in the vicinity of these galaxies, support the hypothesis 
that either a dynamical interaction or a merger with exchange of 
material may have suppressed star formation. Strong AGN/SN 
feedback might also contribute to star-f ormation quenching 
( Springel et aL 20051 Hopkins et al.1 120071) . depending on the 
stellar mass (IKavirai et al. 1120071) . 

This galaxy population was first di scovered, and has been 
primarily studied since, inside clusters (iDressler & Gunn[ ll983: 
Couch & Sharplesll 1987b iTran et al.ll2"003l: iPoggianti et all 1999. 
2009), even if post-starburst galaxies have also been found 
in the field. Numerous studies have been undertaken to study 
the properties of field post-starburst galaxies at low red- 
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Poggianti et al.ll2009h. These studies commonly concluded that 



post-starburst galaxies are a heterogeneous parent population. 

Post-starburst galaxies appear to be rare at any epoch and 
discrepant results about their nature are found in the literature. 
The identifiable sources of discordance are numerous. Besides 
differences in selection criteria and spectral resolution, the main 
source of discrepancy may be the different signal-to-noise ra- 
tio achieved at different redshifts. This implies that we should 
also investigate the dependency of the underlying properties 
on the host galaxies (e.g., the sampling of only the most mas- 
sive or brightest systems). Furthermore, because of the obser- 



vational constraints, the majority of studies cover small areas 
and might probe specific types of environments (large-scale 
structures in pencil surveys, or studies in clusters). In addition 
to the observational limitations, finding post-starburst galaxies 
is also challenging because of t he short duration of the post- 
starburst phase (1— 1.5 Gyr) (see ICouch & Sharpies! d 19871) and 
Barg er et al] (H996)). Because they represent a small fraction of 
the total galaxy population, the systematic effects caused by the 
aforementioned problems require particular attention when in- 
terpreting the results in a general framework. 

In this paper, we present for the first time a mass- 
selected sample of post-starburst galaxies at intermediate red- 
shift in a wide variety of different environments. This mass- 
selected sample and its parent galaxies are extracted from 
the largest multiwavelength survey exist ing up to now, the 
Cosmic Evolution Surv ey, or COSMOS (TScoyille et alJ 120071) 
imaged with HST/ACS dKoekemoer et al.ll2007l) . We provide an 
overview of the main physical properties of these galaxies using 
data of unprecedented quality and completeness in this study. 
We present a first attempt to quantify their cosmic evolution. A 
more detailed study of this subject and its relation to the general 
picture of galaxy evolution over the past 8 Gyr will be presented 
on completion of the survey. 

The present work is organized as follows: the data and sam- 
ple selection are presented in Sect. 2; results are described in 
Sect. 3, and a summary is given in Sect. 4. Throughout this 
work, we assume a standard cosmological model with Q.m = 0.3, 
Q.A = 0.7, and Hq = 70kms~ 1 Mpc~'. Magnitudes are given in 
the AB system. 



2. Sample selection 

2. 1 . The spectroscopic survey: zCOSMOS 

The zCOSMOS project dLillv et al.ll2007h is a redshift survey of 
galaxies in the COSMOS field consisting of two observing pro- 
grams. The "deep part" focuses on high-z galaxies (1.4 < z < 
3.0) using a combined flux-limited (Bab < 25.25) and colour 
selection (these observations are not used in the present study). 
The "bright part" of zCOSMOS project aims to acquire spec- 
tra of about 20,000 galaxies (Iab < 22.5) with a rather uniform 
sampling rate (about 60-70%) across the 1.7 deg 2 COSMOS 
field, a superb spectroscopic redshift success rate (> 97% at 
0.5 < z < 0.8), and high velocity accuracy (~ lOOkms -1 ). 
The bright part is observed with the red R ~ 600 VIMOS MR 
grism covering the spectral range 5500-9500A (2.55A disper- 
sion), which provides an ideal configuration for detecting key 
spectral features in galaxies out to z ~ 1 .2. This paper deals with 
the first "10k" galaxy data set completed to date for the bright 
part of the zCOSMOS survey, which we refer to hereafter as the 
zCOSMOS survey. All details of data acquisition, red uction, red- 
shift m eas urements, and thei r quality can be found in lLillv et alJ 
d2007l) and lLillv et all d2009h . which also describes other impor- 
tant aspects of this survey. 

2.2. The spectral measurements and sample selection 

We measure the spectral features in our 10k sample with an au- 
tomated software package (platefitJVTMOS). H ere we provide a 
brief o utline of the fitting procedure. We refe r to Lamareille et al.l 
(200 6]) for full d e tails, and both iTremonti et all d2004l) and 
Brinchmann et al ] d2004h for further details of the original 
pipeline developed for the higher resolution SDSS spectra. The 
best-fit stellar continuum derived from a grid of stellar pop- 
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Fig. 1. (Top) Distribution of equivalent widths of H6 vs. EW[OlI] 
for all galaxies in the lOk-parent sample in the redshift inter- 
val 0.48 < z < 1.2. This plot shows the criteria adopted to 
select the post-starburst candidate galaxies, i.e., the region be- 
tween H6 A > +3A and EW[OlI]> -3A marked with dashed 
lines. Isophotal contours are in steps of 10% with the faintest 
isophote level starting at the 10% level. (Bottom) Template spec- 
trum of a typical galaxy in a post-starburst phase obtained by 
stacking 69 zCOSMOS post-starburst galaxies more massive 
than \og(MI Mq)> 10.0 (corrected by selection effects as ex- 
plained in Sect. 2.5). 



ulation synthesis models (Bruzual & Chariot 2003, hereafter 
BC03) was subtracted from the observed spectrum. Any remain- 
ing residual was removed by fitting a low-order polynomial to 
the continuum-subtracted spectrum, and all emission lines were 
then fitted simultaneously with a Gaussian profile. Finally, ab- 
sorption features and spectral breaks were measured after sub- 
tracting emission lines from the original spectrum. 



In this paper, we use measurements of the spect ral in- 
dex H6a as defin ed by IWorthev & Ottavianil (1 19971) and 
iBalogh et af] (1 19991) . the equivalent width of the [OlI]/l3727 
doublet (EW[OH]), the equivalent width of the [OIII]/15007 
(EWfOml), and the amplitude of the 4000 A break (D n 4000) 
(Baloghetal. 1999). A careful treatment of spectral fitting 
is mandatory because both emission and absorption lines are 
present in the Balmer lines. The index H6a and the equivalent 
width of the lines are defined such that positive values refer to 
absorption features. 

Authors usually adopt H6a and EW[OlI] features in their cri- 
teria for selecting post-starburst galaxies, and/or Ha depending 
on the sp ectral coverage, a lthough some exceptions exist. For 
example, iTran et all d2004l) used a combina tion of HS and Hy 
indices, in addition to rOHl. lYan et al.l (120061) added an empiri- 
cal limit to EW(H/3) in constraining EW[OlI], and some others 
used, in addition to a threshold on emission lines, the relative 
abundance of A-typ e to K-type stars, i.e., the ratio A/K (e.g., 
lOuintero et al.ll2004l) . In this paper, we adopt a conventional def- 
inition using both the [Oil] emission line and H6a, to be able to 
compare with the majority of other works, although any compar- 
ison is difficult and only possible qualitatively. In a future paper, 
based on the 20k galaxy sample, we will examine the effects of 
different selec tion criteria on a large r sample of H<5-enhanced 
galaxies (e.g. , Le Borgne et al. 2006) and the possible contribu- 
tion from obscured AGNs jYan et al.l2006l) . In the present paper, 
our post-starburst candidate criteria are that the equivalent width 
of HSa is larger than +3A and the equivalent width of the emis- 
sion line [Oil] is larger than — 3A. This is illustrated in the top 
panel of Fig.[TJ where the equivalent width of HSa is plotted ver- 
sus EW[OlI] for the lOk-parent galaxies in the redshift interval 
0.48 < z < 1.2 (the redshift bin explored in this work, which is 
imposed by our instrumental set-up). This plot shows the spec- 
troscopic criteria adopted to select the post-starburst candidate 
galaxies. As is evident in Fig.[T] the region corresponding to our 
k + a candidates is not completely separate from that of the 10k- 
parent populat ion as i t is in the SDSS data (cf. the left panel 
of Fig. 2 by lYan et al.l (120091) 1. This is mainly due to the errors 
associated with the spectral measurements. To evaluate their im- 
pact on the classification scheme, we tested different empirical 
thresholds and applied a confidence level to the spectral fea- 
tures when defining our samples. This approach aims to exclude 
(include) residual signatures of ongoing star formation and to 
test the reliability of our spectral measurements. While the num- 
bers of objects change depending on the thresholds adopted, of 
course, we obtained an optimal compromise between the abil- 
ity to complete a robust, statistical analysis and reliable results 
by including only galaxies with EW[OlI] larger than -3A and 
EW(H(5) larger than 3A at the 2<x confidence level. 

To maintain strong control on occasional poor wavelength 
calibration that can affect in particular the bluest part of our 
spectra, we visually inspected each spectrum of potential can- 
didates. For spectra with doubtful measurements of EW[OlI], 
we imposed an additional limit on the equivalent width of the 
EWfOlII] line. We used a subsample of high signal-to-noise 10k 
spectra and high confidence spectral measurements in which 
both EW[OlI] and EW[OlII] lines were detectable in the spectra, 
to define this limit on the [OIII] lines. From the best-fit model 
solution, we inferred a threshold for the EW[OlII] of — 6.9A. 
Furthermore, from all of our catalogues, we removed both stars 
and broad-line AGNs and adopted a high-quality flag scheme 
for the spe ctroscopic r edshif ts (> 95.5% confidenc e level) (for 
details see iLillv et aD d2007l) and iLillv et all (120091) 1 Although 
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we inspected each spectrum individually, our confidence level in 
the redshift measurements is 100% for our post-starburst sample. 
We also performed quality checks of our post-starburst candi- 
date spectra to verify the impact of residual sky lines on spectral 
features. We examined the reduced VIMOS spectrum for the im- 
pact of fringing effects, which may be significant at wavelengths 
above A > 8,500A be cause of the thinned , back-illuminated 
CCDs used in VIMOS dScodeggio et alj|2005l) . 

The final sample presented in Sect. [3] contains spectra that 
passed all of the aforementioned tests. 

Using the constraints described above, we obtain a sample of 
74 post-starburst galaxies in the flux-limited zCOSMOS bright 
sample in the redshift interval z e [0.48 -1.2]. Three of them are 
at redshifts higher than z > 1 . For one object, the computation of 
a precise absolute magnitude is impossible because of the lack 
of sufficient photometric coverage, and it is excluded from our 
statistical analysis (for details see in Zucca et al. (2009)). 

2.3. Post-starburst spectral template 

Given the difficulties in selectingthis rare population of galaxies, 
we provide a template spectrumj to the community of the post- 
starburst galaxy population obtained by stacking 69 zCOSMOS 
post-starburst galaxies more massive of log(A1/A1o)> 10.0 (cor- 
rected by selection effects as explained in Sect. 2.5). A descrip- 
tion of the steps performed to obtain this template spectrum is 
given in the following text. 

The bottom panel of Fig.Q] shows the template spectrum, 
which exhibits the distinctive spectral features of galaxies that 
have recently experienced a starburst before its star-formation 
was quenched prior to observations. These features can be sum- 
marized as follows: 

i) The strength of the [Oil] nebular emission line, which is 
emitted in HII regions around O and B stars, is almost absent 
(EW[OlI]< -2.5 A) in our template spectrum. Given that [Oil] is 
an indicator of recent star formation (with lifetime of 10 7 yr), this 
property indicates that the star formation during the epoch of the 
observation has stopped. The Ha emission line is another robust 
indicator of ongoing star formation, but cannot be detected with 
our adopted R ~ 600 VIMOS MR-grism observational setup, 
which covers the spectral range 5500A < A < 9500A. 

ii) Another typical spectral feature of this galaxy population is 
the strong absorption features of Balmer lines. Although Ha and 
HyS are often filled in by nebular emission that can make their 
interpretation difficult, HS can be used to estimate the stellar 
ages. From late-type stars to earlier-type stars, the H6 strength- 
ens. This can be achieved if A-stars dominate the main-sequence 
contribution near the epoc h of the obse rvation. 

iii) As noted originally by iRosd dl984l) . the ratio of the equiv- 
alent widths of the Call H/t3968+He to Call KA3934A spectral 
lines is constant in stars later than about F02, but increases dra- 
matically for earlier type stars as the Call lines weaken and He 
strengthens. This spectral property is also illustrated in our tem- 
plate spectrum, where the Call H/13968+He is deeper than the 
Call K/I3934A absorption lines. 

iv) We note the lack of emission from [OIII] and H/3 lines in our 
template spectrum, which indicates that narrow-line Seyferts do 
not represent an imp ortant compone nt of our sample, in contrast 
to other works (e.g. jYan et al .1120091) . 



1 This template can be accessed from the web 
http : / /www . bo . astro . it/~ daniela/ template .fits 
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Fig. 2. Stellar mass and B-band absolute magnitude as a func- 
tion of redshift of the zCOSMOS bright sample (points, gray 
symbols). The (violet) triangles show the location of post- 
starburst galaxies. The light and dark gray areas show the selec- 
tion in stellar mass adopted in sample MCOR and sample MST, 
respectively. The dashed (black) line in the upper panel illus- 
trates the completeness in mass-to-light ratio for the global pop- 
ulation (for comparison, it is plotted with a dotted, red line the 
completeness limit for the early-type population). The dashed 
(black) line in the bottom panel is plotted as reference to il- 
lustrate that our galaxies of sample MST are brighter than the 
evolving-luminosity limit at Ms < -20.1 - z- 



The combined use of the two SFR estimators (H6 and [Oil]), 
which are sensitive to stars of different stellar masses and life- 
times, allows us to obtain in the next section a quite precise cos- 
mic clock of their recent (< 1 Gyr) star-formation history, and to 
characterize their physical properties. 

2.4. Stellar masses 

Since we studied galaxies that are understood to have expe- 
rienced a recent secondary burst of star formation, we used 
Chariot & Bruzual models with the inclusion of secondary bursts 
(Chariot & Bruzual in prep., or CB07). These models have been 
frequently used to model the complex star-formation history 



lly man a p 
aLl l2003at ll 



Brinchm ann et all [2004; 



of a galaxy more re alistically than a purely exp onential ly de- 
clining law ( Kauffmann et 
ISalimetal.1120051) . 

Stellar masses were derived by fitting stellar population syn- 
thesis models from the library of CB07 to the broad-band optical 
(CFHT: u, i, Ks; Subaru: B, V, g, r, i, z; Capak et al. 2007) and 
near-infrared (Spitzer/IRAC: 3.6/im, 4.5pm; Sanders et al. 2007) 
photometry by minimizing the chi-square test statistic for each 
galaxy. The estimate of the stellar mass provided by the best-fit 
model spectral energy distribution is the current mass contained 
in stars in each galaxy at the epoch of the observation obtained 
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by integrating the star formation history over the galaxy age and 
subtracting the cumulative mass lost during stellar evolution. 

The measurement of stellar mass included some assump- 
tions. We used a Chabrier initial mass function (Chabrier 2003) 
with a lower and upper cutoff at 0.1 and lOQyVto, a d ust attenua- 
tion described by the model of lCharlot & Falll d2000l) . and metal- 
licities between 0.1 and 2 Z©. The star-formation history was 
modelled by two components. The first component was a contin- 
uous star-formation model with an exponentially declining law 
of the form SFR(t) oc exp (-t/r), where timescale t and age t are 
in the range r/[Gyr] =[l,oo] and t/[Gyr] = [0.1, 20]. The second 
component was represented by random bursts superimposed on 
the first component that had an equal probability of occurring 
throughout the life of the galaxy. 

Of course the results depend on the adopted assumptions and 
models used. For example, the assumption of a Chabrier initial 
mass function produces a consistent difference in the values of 
the stellar mass at fixed age that is within a factor of 1 .7, but with 
a small di spersion, of those d erived with the Salpeter IMF pre- 
scription (Poz zetti et al.ll2.Q07h . The inherent uncertainty in the 
initial mass function adopted is not a global limitation for re- 
sults based on stellar masses given the constancy of this rigid 
shift over a wide range of star formation histories, although to 
limit the various degeneracies a large set of grids of stellar pop- 
ulation synthesis models were used. Further details of the s tellar 
mass measurements can be found in lBolzonella et al] d2009l) and 
IPozzetti et all (l2007ll2009t) . 

2.5. Post-starburst representative sample 

Since our original selection is based on a flux-limited survey, we 
used subsamples of galaxies complete in terms of stellar mass 
(adopting two independent approaches) to avoid selection bi- 
ases. 

In the first approach, we adopted a conservative limit to 
the stellar mass of \og(M/M B )> 10.8 in the redshift bin z e 
[0.48 - 1.0] (from now on, this stellar Mass STringent crite- 
rion refers to sample AiST). This stringent mass limit was de- 
fined using Mock catalo gues from the Millennium simulation 
dDe Lucia & Blaizodl2007l) of the zCOSMOS bright sample (see 
Meneux et al. (2009J) for an exhaustive explanation on the sub- 
ject). We selected 18 post-starburst galaxies for this conservative 
mass-complete sample, or sample AiST. 

To maximize the quality of the statistical analysis, we also 
constructed a second sample, or sample MCOR, including 
galaxies more massive than log(A1/Al0)> 10.0 in the redshift 
bin z € [0.48 - 1 .2], and accounting for under-represented galax- 
ies at high z by correcting with statistical weights and a V tot /V max 
formalism ( Schmidt 1968). The total volume sampled by the sur- 
vey within the redshift bin explored is denoted by V tot . The vol- 
ume V max is the maximum volume in which each galaxy of a 
given Iab magnitude is still observable. We applied the statistical 
weights adopted in the studie s of luminosity and mass functions 
of the zCOSMOS sample dZucca et all 120091; iBolzonella eTal] 
2009: IPozzetti et al. 2009]). These weights are the product of the 
inverse of the target sampling rate (TSR) and the inverse of the 
spectroscopic sampling rate (SSR). The target sampling rate is 
the ratio of the objects spectroscopically observed to the total 
number of objects in the parent photometric catalogue. The spec- 
troscopic sampling rate is the ratio of the number of observed 
objects with reliable redshifts to the total number of observed 
objects, calculated in bins of apparent magnitude. The apparent 
magnitude dependence takes into account the decline towards 
fainter magnitudes, in our ability to measure a redshift. A more 



complex scheme that introduces a redshift dependence of SSR 
does not alter the final results appreciably in the redshift bin ex- 
plored here (see Bolzonella et al. 2009). In the sample MCOR, 
the mass completeness limit, which depends on both redshift and 
mass-to-li ght ratio, also takes account of the colou r-magnitude 
relation dPozzetti et al.|[2009i: IBolzonella et al.ll2009l) . 

The second approach permits us to increase the number of 
objects (from 18 to 69) in our analysis and explore a wider red- 
shift bin (z e [0.48 - 1 .2]) than that of z e [0.48 - 1 .0] for sample 
AiST. However, we always compare the results obtained with 
this MCOR sample with those obtained with the complete MST 
sample to validate further our conclusions. 

Figure|2]shows the stellar mass of the zCOSMOS bright sam- 
ple as a function of redshift and, overplotted, the data of the post- 
starburst galaxies in our mass-complete samples (with violet, tri- 
angle symbols). We note that the conservative limit imposed by 
Mock catalogues agrees well with the completeness in the ob- 
served mass-to-light ratio for the global population. This limit is 
similar to the more conservative limit obtained for the early-type 
population (defined using the spectrophotometric classification 
described in Zucca et al. 2009). The statistical weights and mass 
completeness adopted are the values u s ed in t h e luminosity and 
mass functi on studies by IZucca et aTl (2009), Bolzonella et alj 
(2009), and IPozzetti et alJ d2009l) ."ln the same figure, we also 
plot the B-band absolute magnitude as a function of redshift. It 
is important to note that our sample AiST, which is complete in 
terms of stellar mass such that log(Al/A1o)> 10.8, is brighter 
than the evolving-luminosity limit at Mg < -20.1 -z (illustrated 
with the dashed, black line overplotted in the left panel of Fig.[2j». 

2.6. Comparison samples 

We compiled two comparison catalogues of strongly star- 
forming and quiescent galaxies. For these galaxies, we also 
adopted consistently the same stellar mass completeness crite- 
ria, redshift flags, and confidence levels on the spectral measure- 
ments as used in defining the class of post-starburst galaxies. 

The star-forming and quiescent galaxy classes were selected 
using a pure spectral classification. We combined measurements 
of the amplitude of the 4000 A break, taken as an age esti- 
mator of the underlying stellar population, and the equivalent 
width of the [OlI]/l3727 line, assumed to be an indicator of 
the star-formation activity. Thus, quiescent galaxies were de- 
fined to be galaxies with a zero, or negligible, star-formation 
rate inferred using the EW[OlI], i.e., EW[OlI]> -3A, and an 
old underlying stellar populat ion described by their lar ge 4000A 
break (D n 4000> 1.5) (e.g.. iKauffmann et all d2003bh l. We de- 
fined a sample of strong EW[OlI] emitters by constructing a 
catalogue of objects with EW[OlI]< -15A and small 4000 A 
break (D n 4000< 1.5). From the 10k sample, we also compiled 
the lOk-parent catalogue of galaxies following the selection in 
stellar mass adopted for sample AiCOR and sample AiST, and 
all other criteria listed above (from now on when we use the ter- 
minology "all objects", we refer to the lOk-parent catalogue). 

Our selected samples are displayed in Fig. [3] where the selec- 
tion criteria used to classify galaxies in our comparative analysis 
are shown. This figure illustrates how the three classes of galax- 
ies (highly star-forming, quiescent, and post-starburst galax- 
ies) populate the plane of the amplitude of 4000 A break and 
EW[OlI]. Small points and isophotal contours in steps of 15%, 
with the faintest isophote level starting at 15%, show the 10k- 
parent galaxies level. Using criteria adopted for sample AiCOR, 
we obtain a comparison sample of 295 quiescent galaxies and 
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Fig. 3. The amplitude of 4000 A break is plotted as a func- 
tion of EW[OlI] to show the selection criteria used to classify 
the classes of star-forming (square, blue symbols) and quiescent 
(circle, orange symbols) galaxies. Triangle (violet-coded) sym- 
bols show the post-starburst galaxies. Small points (in black) 
show the lOk-parent catalogue (isophotal contours are in steps 
of 15% with the faintest isophote level starting at 15% level). 
The 4000 A break histograms of the lOk-parent sample (dashed, 
gray line) and of the post-starburst galaxies (solid, violet line) 
are also shown along the x-axis. The vertical lines indicate the 
location of the cut in D n 4000 (at 1.5). The horizontal lines indi- 
cate the location of the two extreme thresholds at -3A and -15A 
in EW[OlI]. Arrows are upper limits to EW[OlI] values. Typical 
errors in the spectral measurements are plotted as errorbars at 
the bottom of the panel. 

174 star-forming galaxies, while we have 178 quiescent galax- 
ies and 45 star-forming galaxies in sample MST. The lOk-parent 
galaxies number 1508 following the criteria of sample MCOR, 
and 665 for sample MST. 

3. Results 

We present the main physical properties of the mass-selected 
post-starburst galaxies identified in the zCOSMOS survey in at- 
tempting to obtain a coherent picture of their origin in the con- 
text of galaxy formation. Our analysis refers to galaxies of the 
statistically corrected, mass-complete sample at z e [0.48 - 1 .2] 
(sample MCOR) and the conservative mass-limited sample at 
Z e [0.48 - 1 .0] (sample MST). 

3.1. Robustness of our sample: dusty starburst hypothesis? 

Galaxies with spectral characteristics similar to those used to 
define a post-starburst galaxy are often interpreted to be a sys- 
tem in which star formation has suddenly and almost com- 
pletely stopped within the past few Gyrs. This cessation of star- 
formation is assumed to explain the lack of emission lines in the 
galaxy spectrum. 

An alternative explanation, however, is that these galaxies 
do have emission, but that this is extinguished by a significant 



amount of dust. Galaxies with these characteristics should be 
dusty, starburst galaxies (and misclassified post-starburst ob- 
jects). This is a hypothesis that should be verified especially 
when dealing with galaxies understood to have experienced a 
recent, strong burst of star formation. OB stars responsible for 
the [Oil] emission, as other younger stellar generations, are 
more strongly obscured by dust than older stellar populations. 
This age selective dust extinction enables the true strength of 
the Balmer absorption lines to be measured, but only a frac- 
tion of the emission to be de tected in stronger starburst galaxies 
(Poggianti et al.ll2~00lL 12009). To ascertain whether our sample 
of post-starburst galaxies may be contaminated by dusty star- 
burst galaxies, we performed different tests. First, we explored 
the dusty starburst scenario by computing both the extinction 
values and star-formation rate by SED fitting. Secondly, taking 
advantage of the available multiwavelength obs ervations, we ex- 
plored the X-ray and 1 .4 GHz radio data sets dSchinnerer et al] 
120071: ICappelluti et alj|2009l:lBrusa et alj|2009l) . 

1) Dust extinction and SED fitting SFR. We found that post- 
starburst galaxies have in general low values of extinction as esti- 
mated by fitting CB07 model SEDs (median value < Ty >- 0.6). 
The median SFR is 0.6A1oyr _1 , and the 25% and 75% quartiles 
of the distribution range from O.OlMoyr 1 to 4.2M yr~ 1 . The 
distributions of both Ty and SFRs for post-starburst galaxies dif- 
fer at a high confidence level from the same distributions for 
star-forming galaxies (< ry >~ 2 and < SFR >~ 2IAI0). While 
post-starburst galaxies exhibit a similar cumulative ry distri- 
bution to quiescent galaxies, the latter galaxies have an even 
lower rate of star-formation activity, as expected when using 
the used spectroscopic criteria described in Sect. 2 to identify 
them (< r v >~ 0.4 and < SFR >~ 0.07 M Q ). For only one 
post-starburst galaxy do we infer from SED fitting a high rate 
of star-formation activity (~ 140M e yr~ l ) and a high extinction 
(ry ~ 2). This object is one of the two post-starburst galaxies de- 
tected at both X-ray and radio 1 .4 GHz frequencies, even if the 
radio detection is only at the 3<x confidence level (see below). 

2) 1.4 GHz radio frequency. We investigated the dusty star- 
burst hypothesis further by cross-corr elating our sample with 
the 1.4GHz VLA -COSMOS survey dSchinnerer et all 120071: 
iBondi et alJ2008l) . If it is not produced by AGN emission, the ra- 
dio flux is indicative of the current star-formatio n activity (over 
short timescales of 10 8 yr) unbiased by dust (see ICondonl f 1 992, 
for a review). 

We found that one post-starburst galaxy has a counterpart 
at 1.4 GHz radio frequency with a detection confidence level 
of 7.5<r and five others have possible 1.4 GHz emission at 
3<x confidence level. For the entire sample of 74 post-starburst 
galaxies, we therefore have a detection rate of 1.3±1.3%, if we 
adopt the threshold o f 4cr for assuming a ra dio detection to 
be tru e (as adopted in ISchinnerer et al] d2007l) and IBondi etall 
(2008)). This detection rate increases to 8.1+2.7%, if we adopt 
the less conservative threshold of 3cr for radio detection. 

To investigate the mean radio properties of the post-starburst 
galaxies, we performed a stacking analysis of the 1.4 GHz radio 
images. From the stacking analysis, we excluded the six post- 
starburst galaxies with a possible radio counterpart plus three ad- 
ditional sources that exhibit, within 2.5" of their optical position, 
either radio emission more significant than 3<x that is not asso- 
ciated with the post-starburst galaxy, or a negative radio peak of 
higher significance than 3<x. For the remaining 65 post-starburst 
galaxies, images of 40 x 40 pixels (corresponding to 14" x 14") 
were extracted around the optical position. Finally, all the 65 im- 
ages were combined to obtain a median image with a lcr radio 



D. Vergani et al.: Post-starburst galaxies in COSMOS field 



7 



noise of *2 //Jy. A marginal detection of « 3<t confidence level 
was found at a distance of 1.4" from the centre of the stacked 
image with a radio flux of ^7 yuJy. To access the reliability of 
this detection, we repeated the stacking analysis for 5 different 
samples, each of them consisting of 75 random positions. The 
criterion used to exclude from the stacking, post-starburst galax- 
ies with a radio detection (or a negative radio peak, see above) 
was applied similarly to the simulated samples. All 5 stacked 
images from the simulated sample had a lcr radio noise in the 
range 1.4-1.8 /Jy. No radio detection was identified for any of 
the simulated stacked images. 

Therefore, assuming the marginal detection in the stacked 
image to be real and assuming its flux of 7 yuJy to be the mean 
radio flux of radio-quiet post-starburst galaxies, we can estimate 
their mean star-formation rate, or more realistically its upper 
limit since some objects may be AGNs. We obtaine d an upper 
limit of 8.2A1 Q yr _1 by adopting the calibration by iBell et al.l 
d2003h . The radio luminosity Li 4ghz was calculated by assum- 
ing a median redshift of Zmedian = 0.7014. 

We then explored the physical properties of the six (> 
3cr) radio-detected post-starburst galaxies. Radio-detected post- 
starburst galaxies are slightly more massive and brighter than 
non-radio-detected post-starburst galaxies (logvVt/Alo! 4gh z = 
11.22 ± 0.09 versus log Al/Me no _, 4GH z = H-08 ± 0.05; 
M B ,i.4GHz = -22.64 + 0.37 versus M B ,„ -i.4GHz = -22.44 + 0.14). 
Post-starburst galaxies with a counterpart at 1 .4 GHz also show 
slightly bluer colours ((U-B)i.4ghz = 1-10 + 0.1) than non-radio- 
detected post-starburst galaxies ((U-B) no _i 4ghz = 1.18 + 0.03). 
Although these differences are not statistically significant per 
se, their properties and low detection rate suggest that radio- 
detected post-starburst galaxies might represent the initial phase 
in the evolution of post-starburst galaxies, close to the peak of 
the starburst (the timescale of radio emission is an order of mag- 
nitude shorter than that of the k + a phase). 
3) X-ray counterparts. To characterize the multiwavelength 
properties of post-starburst galaxies, we also searched for their 
possible X-ray counterparts using the available data provided 
by XMM-Newton for the entire 2 deg 2 of the COSMOS field 
(Cappelluti et al. 20 09; Brusa et al. 2 009) and Chandra over 
the inner 0.9 deg 2 dElvis et al.1 [2009) at considerably deeper 
X-ray flux levels. While only one post-starburst galaxy is de- 
tected by XMM-Newton for a matching radius of 6", five X-ray 
sources (four "new" plus the XMM-Newton one) are detected by 
Chandra (within 4" because of the sharper PSF) at a detection 
likelihood >10.8 in at least one X-ray band (see Elvis et al. 2009 
and Puccetti et al. 2009 for further details) among the 54 post- 
starburst galaxies located in the field of Chandra mosaic. Visual 
inspection of all of these sources supports their likely physi- 
cal association with the post-starburst galaxies under investiga- 
tion. All of the five post-starburst galaxies detected by Chandra 
are in regions of moderate-to-high exposure time (>100 ks), 
where more than one Chandra pointing is available (Elvis et al. 
2009). However, the limited number of counts 9^4-0) in the 
0.5-7 keV band prevents us from performing an adequate X-ray 
spectral analysis. X-ray counts were converted into fluxes (us- 
ing a power-law model with photon index F = 1 .4, as adopted 
in the Chandra-COSMOS catalog; Elvis et al. 2009), and then 
into 0.5-10 keV rest-frame luminosities (given the source spec- 
troscopic redshifts of 0.67-1.10) of l.l-8.6xl0 42 erg s -1 . 

At least for the three sources with the highest X-ray counts 
in the hard band (above 2 keV), this luminosity range suggests a 
significant contribution from AGN emission. The assumption of 
a photon index of F = 1 .9, which is typical of unobscured AGN 
emission, would decrease X-ray fluxes and luminosities by 25%. 



The optical spectra of these three sources with the highest X-ray 
counts in the hard band, might indicate that they are AGN host 
galaxies. 

Of these five X-ray sources, two ha ve a clear as sociation with 
sub-mJy sources in the VLA catalog dBondi et al. 2008D, one of 
which has the highest SFR of all sources and dust extinction 
inferred from SED fitting. 

To place constraints on the properties of the X-ray unde- 
tected post-starburst galaxies, X-ray stacking analysis was ap- 
plied to the remaining 49 sources, extracting source counts from 
circles of 5" radius and background counts from a region at 
least 10" away from the source position. X-ray sources present 
in the Chandra source catalog (Elvis et al. 2009) were care- 
fully masked. To limit the problems with the PSF at large off- 
axis angles, we extracted source counts using only observations 
where t he source position w as l ocated within 8' o f the field cen- 
tre ( see iMivaii et all d2008l) and lKim etail d2009l) for a descrip- 
tion of the tool adopted here). We did not find any significant 
signal from the stacked post-starburst galaxies using Chandra; 
the 3cr flux upper limit to the stacked 0.5-8 keV counts is 
« 2 x 10~ 16 cgs, corresponding to a rest-frame 0.5-8 keV lumi- 
nosity of « 4.9 x 10 41 erg s _1 at the average redshift (( z) ~ 0.73) 
of the stacked sources. This upper limit to the stacked signal is 
consistent with emission from normal and starburst galaxies. 



iDewangan et al.l (|2000) reported a detection of an active nu- 
cleus in one post-starburst galaxy selected during a search for ul- 
trasoft sources in a ROSAT campaign, but no radio emission was 
detected within a radius of 15" in the NRAO/VLA and FIRST 
sky surveys. In previous studies, simila r counter-checks led to 
a detection rate of the order o f 1-15%. iMiller & Owed ([2002 ) 



19961) . detected 



using the galaxies selected by IZabludoff et al.l 
2 in 15 post-starburst systems in radio continuum (correspond- 
ing to a star-formation rate of a bout 2 - 6 yVlp yr -1 ), and 1 in 
56 objects were detected by Bla ke et al.l d2004l) (at 7 A1 yr _1 ). 
Using C handra observations of th e Extended Growth Strip from 
DEEP2, Georgakak is et alJ d2008) found an association between 
post-starburst galaxies at z ~ 0.8 and AGN, supporting a mod- 
erate levels of AGN obscuration. Although these results are dif- 
ficult to compare because of different selection criteria and dif- 
ferent mass and luminosity ranges, these results have been used 
to claim that at least some "so-called post-starburst" galaxies are 
still actively forming stars, implying that a selection based of 
the blue-part of the optical spectra is imprecise. Some authors 
have shown that optical spectra typical of a post-starburst galaxy 
in th e blue part of their spec trum, have substantial Ha emission 
(e.g., Liu & Kennicutt 1995). We fully agree with the statement 
by e.g., IMiller & Owed (120021) that the minimum requirement 
for defining a system in their k + a, or post-starburst phase is the 
coverage of both blue- and red-parts of the optical spectrum, or 
complementing the optical data with other wavelength data sets. 



Based on our complementary radio-continuum and deep X- 
ray data, we found that a negligible fraction of our post-starburst 
galaxies, selected by using blue diagnostics of the optical spec- 
tra, have possible residual star-formation activity, and a small 
percentage shows a contribution from AGN emission. Thus, the 
main conclusion of our present study is that we have selected a 
robust sample of post-starburst galaxies. 
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Fig. 4. Stellar mass, B-band absolute magnitude, and stellar mass-to-light ratio histograms are shown for each class of galaxies. The 
cumulative distributions of each sample are also plotted. Median values of post-starburst galaxies and lOk-parent sample are plotted 
as solid and dashed lines, respectively. Blue and red arrows are median values of star-forming and quiescent galaxies. The bars 
assigned to each median value correspond to Icr error computed with resampling technique. Galaxies classified as post-starburst 
(violet), quiescent (red), star-forming (blue), and the lOk-parent (black) sample follow the definition described in Sect. 2.5. 



3.2. Basic properties: luminosities, stellar masses, and 
colours 

The B-band absolute magnitude distribution plotted in Fig. [4] 
shows that post-starburst galaxies of the conservative sam- 
ple MST populate the brightest tail of our lOk-parent galaxy 
population. The post-starburst galaxies more massive than 
log(A1/Alo)> 10.8 are brighter by ~ 0.7 - 0.8 mag than the 
galaxies in both our comparison classes. A similar analysis of the 
sample MCOR galaxies confirms that the post-starburst galaxy 
sample is systematically brighter by ~ 0.4-0.5 mag than the me- 
dian brightnesses compared to median values of the other galax- 
ies. Both these results are statistically significant at more than 
the 3cr level. One of the reasons for this B-band brightness ex- 
cess is probably a recent burst of star formation in post-starburst 
galaxy candidates, as commonly concluded in past studies. 

Post-starburst galaxies also populate the massive tail of the 
lOk-parent galaxies (see Fig.|4]i. Their mass distribution differs 
from that of the lOk-parent galaxies at ~ 91% confide nce level 
using the Kolmogorov-Smirnov test dPress etalJll992h . and our 
post-starburst candidates are as massive as quiescent galaxies. 
Our different samples selected using the MCOR criteria show 
the typical increasing trend in stellar mass and mass-to-light ra- 
tio expected on the basis of their galaxy classification from star- 
forming galaxy population to quiescent galaxies. Given the sim- 
ilar magnitudes and stellar masses spanned by the star-forming 
and quiescent galaxies selected using the MST criteria, this trend 
is not visible in the corresponding panels of Fig. [4] Post-starburst 
galaxies appear to have a brighter B-band luminosity because of 



a combination of a high stellar mass and relatively low mass-to- 
light ratio. Among the most massive galaxies (criterion MST), 
there is a trend for post-starburst galaxies to be the most mas- 
sive objects, more massive than quiescent galaxies (at Icr) and 
star-forming galaxies (at ~ 2(f). Since our post-starburst galax- 
ies are already so massive, we may conclude that their progen- 
itors cannot be common star-forming galaxies that have assem- 
bled their stellar mass by means of a continuous, smooth star- 
formation history and that have stopped, for whatever reason, 
forming stars. Indeed to justify the observed stellar mass as- 
sembly of a post-starburst galaxy with a typical stellar mass of 
log(M/M )= 10.58 (see Fig.g}, a specific-SFR (SFR/M) of 
~ 1 .2 x 10~ 9 (or SFR~ 20 Af© yr _1 ) is required for a galaxy with 
a stellar mass of log(A1/A40)= 10.23 (where 10.23 is the median 
stellar mass found for our star-forming galaxies). This specific- 
SFR value is typical of an object with a high star-formation ef- 
ficiency. Furthermore, we can also exclude the dry merging hy- 
pothesis, i.e., the merging of galaxy progenitors in the absence 
of gas, as the strength of the H6 absorption lines implies a recent 
burst of new stars, i.e., that fresh gas should have been recently 
available. Therefore, the possible progenitors might be the al- 
ready massive, less numerous star-forming galaxies at each cos- 
mic time, or/and the most actively star-forming sources that in a 
short period of cosmic time can assembly sufficient stellar mass: 
in other words, the k + a's progenitors might be galaxies with 
the highest values of specific star formation rates at any epoch, 
assuming that the mass assembly is a process sustained by star 
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formation activity at the observed rate, i.e., no merger events are 
invoked. 

Quiescent and star-forming galaxies populate the red se- 
quence and the blue cloud, respectively, as shown in Fig. [5] for 
sample MCOR. The peak of the distribution for quiescent galax- 
ies is at (U - B) ~ 1.3, and star-forming galaxies are distributed 
around a locus of (U - B) ~ 0.8 in sample MCOR (colours 
are slightly redder for star-forming galaxies in sample MST by 
~ 0.2 mag, due to the higher stellar mass selection). The un- 
derpopulated region in-between the two distributions, i.e., a re- 
gion called the green valley, is the locus within which our post- 
starburst galaxies are located. In the colour-magnitude diagram, 
they are on average (U - B) ~ 0.25 redder than star-forming 
galaxies, and ~ 0.1 bluer than quiescent galaxies. Although the 
median distribution of colours appears to imply that the green 
valley is populated by post-starburst galaxies, they span a wide 
range of colours. This property is confirmed in Fig. [3] where the 
4000 A break histograms of the lOk-parent sample and post- 
starburst galaxies are shown along the x-axis (with a solid, vio- 
let line and dashed, gray line, respectively). The 4000 A break of 
post-starburst galaxies ranges between quite small values (typi- 
cal of star-forming galaxies) and quite high ones (characteristics 
of quiescent, passive galaxies), or 1.2 < D„4000 < 2, with a 
peakatD„4000 - 1.5. 

Based on this result, selecting post-starburst-candidate 
galaxies such as those that exhibit stronger Balmer absorption 
lines than expected for thei r 4000 A break, e. g., in a way similar 
to the criterion adopted bv lWild et alJ d2009l) . should underesti- 
mate their true fraction. Our result might justify the low detec- 
tion of this class of galaxies in the VVDS field. Wild et al. find 
only 5 post-starburst galaxies in the VVDS field, scaling from 
their quoted 18 HS strong galaxies with nebular emission lines 
to 5 properly-defined post-starburst galaxies (i.e., those objects 
among the 18 detected ones without nebular emission lines). 

Our result indicates that post-starburst galaxies, which pop- 
ulate the green valley for a very short time, provide a suitable 
galaxy population to occupy the gap between blue and red galax- 
ies in the colour-magnitude relation observed at any redshift 
from the local Universe up to z ~ 1. — 1.5. Although k + a galax- 
ies may represent some of the galaxies in the green valley, we do 
not however propose that they constitute the entire population of 
green-valley galaxies. 

3.3. Morphologies and structural parameters 

We investigated the morphological properties of post-starburst 
galax ies using the high spatia l resolution HST/ACS F814W im- 
ages (Koe kemoer et alj |2007). Galaxy morphologies contain the 
imprints of stellar populations, gas content, and dynamical struc- 
tures. Thus, if an event such as merging has occurred in the re- 
cent history of a galaxy, its signature should be observed, given 
that the timescale of gravitational disturbances induced by merg- 
ing is similar to that of the post-starburst phase. 

We used the ZEST morphological classification of 
IScarlata et alJ (l2007al) based on five non-parametric diag- 
nostics to define early, late, and irregular morphological classes 
of all studied galaxies in this work. We also explored in more 
detail properties such as the asymmetry index, the degree of 
irregularities, and the concentration of the light. 

The upper panel of Fig. [6] illustrates that a similar frac- 
tion of post-starburst galaxies in the sample MCOR are bulge- 
dominated (E) and disky-dominated (L) (see panel B). How 
the galaxies are distributed between the so-called early- and 
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Fig. 5. Distribution of the (U-B) colour versus the stellar mass 
for post-starburst galaxies (triangle, violet colour-coded sym- 
bols) and comparison samples (quiescent galaxies, circle, red- 
coded symbols, and star-forming galaxies, square, blue-coded 
symbols). Isophotal contours are in steps of 10% with the 
faintest isophote level starting at the 10% level. The vertical solid 
(dashed) line indicates our mass-completeness at log(A1/A1o)> 
10.0 (> 10.8). 



late-type divisions is independent of their H6 strength shown 
in panel A of Fig. [6] There is, instead, an increase with H6 
in the fraction of post-starburst galaxies in irregular systems. 
In panel C, we show the relative fractions of the three dif- 
ferent morphological types (bulge-dominated, disky-dominated, 
and irregulars) that are post-starburst galaxies compared to 
the same fractions for the samples of quiescent and star- 
forming galaxies. As expected, star-forming galaxies are mostly 
disky-dominated with ~ 20 - 30% having irregular morpholo- 
gies, while the majority of the quiescent galaxies are bulge- 
dominated. Morphologically, post-starburst galaxies are a het- 
erogeneous population, although we note that the most massive 
post-starburst galaxies (sample MST) are predominantly bulge- 
dominated systems (see lower panels of Fig.|6j and have a sim- 
ilar distribution among morphological types as quiescent galax- 
ies. It is interesting to note that in the sample MST, there is no 
post-starburst galaxy with an irregular morphology, but we also 
note the low number statistics in the current sample MST. 

In some studies, post-starburst galaxies are found to be 
mostly bulge-dominated and/ or have a morphology that is con- 
sistent with being earlv-tvpe dGoto & SPSS collaborationl 2003 ; 

DtO 



Ouint ero et"aTll2004t iTran et al.ll2004t iBlake etalj 12004 iGotc. 
12005b iBalogh et alJ (20051: iPoggianti et alJ l2009tT and in some 
others have a high i ncidence of disks an d spiral arms (e.g.,[F ranx 
1 1993b [ Couch et alj|1994t iDressler et a lJ I 19941 ICaldwell & Rose! 
1 19971: Dressier et alJll999h . For example. lOuintero et all d2004) 
visually inspected a sample of 160 galaxies from the SDSS at 
0.05 < z < 0.20, and classif ied 90% of them as bulge-dominated 
systems. Yan g et aD (120041) used the term "diverse morpholo- 
gies" to describe five high-reso lution HST images of k +a galax- 
ies selected from the sample of Zablud off et al. (Q996), although 
we note that their selection did favour the bluest systems, i.e., 
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the ones with smallest D n 4000 that have experienced most re- 
cently a significant starburst episode. Their detailed study using 
bulge-to-disk two-dimensional fitting decomposition of four of 
them identified one barred SO, two disky systems, and one un- 
classifiable galaxy. Their inferred central light distributions were 
found to have a profile similar to, but more luminous , than that 
of a ty pical early-type power-law. At higher redshifts. lTran et al.1 
(2004) found more bulge-dominated systems in their k+a galaxy 
sample than in field galaxies. The sample was selected from 
a program designed to study galaxy evolution in clusters, and 
k + a galaxies were selected using a conservative constraint on 
both the H6 and Hy Balmer indices and according to the absence 
of [Oil] emission line . A study based on vis ual classification of 
ACS/HST images bv lPoggianti et al.l <|2009) identified a higher 
fract ion of SO and Sa typ es among field post-starburst galaxies 
than iDressler et al.l dl999h who instead found a higher fraction 
of later types among post-starburst galaxies. 

Based on the morphological results obtained for the present 
post-starburst galaxy sample, we agree with the widely estab- 
lished result that post-starburst galaxies are morphologically a 
heterogeneous population. We also found an increasing frac- 
tion of early-type morphologies at increasing stellar mass in 
agreement with th e well known mass-morphology relation (e.g., 
iTasca et alJl2009T) . We might find in the future with larger data 
sets spanning a wider range of galaxy properties, that the flux 
limits and other selection effects might play a critical role in 
the morphological classification, e.g., favouring an early-type as- 
signment for brighter/massive systems as found in this work. 

Apart from being heterogeneous in their morphologies, the 
second common feature used to describe post-starburst galax- 
ies by most authors is the large fraction of asymmetries and 
dumpiness properties associated with them (e.g., fine structure, 
and tidal tails). We quantify the degree of asymmetry in the 
light distribution using the asymmetry coefficient computed as 
the difference between the image rotated through 180 degrees 
and the original fram e following the usual prescrip t ion (e.g ., 
IConselice etal]|2000h (for details see IScarlata et all d2007bl)\ 
Figure|7] shows that for post-starburst galaxies in the sample 
MCOR, the median and range of the asymmetry index are in- 
termediate between values for quiescent galaxies (and the 10k- 
parent sample) and star-forming galaxies. Comparing the asym- 
metry distributions of post-starburst galaxies with those of star- 
forming and quiescent galaxies, they both differ at higher than 
the 99.9% KS-test confidence level. Similar results are obtained 
for the light concentration, which is defined to be the ratio of the 
radius including 80% to that including 20% of the galaxy light 
and measures the central density of the galaxy light distribution. 
In Fig. |7] we show that the light concentration indices of post- 
starburst galaxies are typically in-between the values for qui- 
escent and star-forming galaxies. These differences in both the 
morphological classification and the parameters of asymmetry 
and light concentration of post-starburst galaxies and that of both 
star-forming and quiescent galaxies, suggest that post-starburst 
galaxies represent an intermediate stage of galaxy evolution. 

Using the ZEST scheme, the distortion of the light distribu- 
tion is quantified by the dumpiness index. This index is avail- 
able for galaxies classified as early- and late-types, i.e., not for 
irregular galaxies, which by definition are very distorted and 
for which the index might be associated with only the max- 
imum degree of dumpiness. The flags for the dumpiness in- 
dex range from zero, indicative of a smooth light distribution, to 
three, corresponding t o a clumpy light distribution (see Table 2, 
IScarlata et alj (l2007al) ). We found that the majority of quiescent 
galaxies have a smooth light distribution (86% have flags of 
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Fig. 6. Morphological classification of post-starburst galaxies for 
sample MCOR (top) and sample MST (bottom). In panel A we 
show the fraction of the three different morphological types in 
post-starburst galaxies (bulge-dominated (E), disky-dominated 
(L) and irregulars (I)) as a function of the strength of H6. Similar 
plot is shown in panel B for post-starburst galaxies (indepen- 
dently by the H6 strength). Panel C shows the fraction of the 
three different morphological types in quiescent galaxies (or- 
ange) compared to the sample of post-starburst (violet) and star- 
forming galaxies (blue), respectively. The error bars represent 
the Poissonian errors. 



or 1), whereas star-forming galaxies tend to have a classification 



D. Vergani et al.: Post-starburst galaxies in COSMOS field 



11 



consistent with having a clumpy structure (70.2% have flags 2 or 
3). Post-starburst galaxies have a frequency distribution of con- 
centrations that overall is similar to that of quiescent galaxies. 
This is indicative of a smooth light distribution and consistent 
with a picture in which the star formation experienced in the re- 
cent past by these galaxies, has ultimately faded to produce a 
smoother morphology. 

Using Hubble Space Telesc ope (HST) WFPC2 im ages, and 
in agreement with our results, Caldw ell et al ] (119991) detected 
smoother morphologies in six post-starburst galaxies than in the 
starburst galaxies of their sample. In contrast, the latter galax- 
ies had centrally concentrated star formation (we note, however, 
weak [Oil] emission in two of the six objects that ide ntify them 
more appropriately with the H<S strong galaxy class). rTran et alj 
(2004) detected both centrally concentrated light emission and 
clumpy regions in the HST WFPC2 images of five post-starburst 
galaxies. Although no conclusive results can be drawn because 
of the widely ranging properties of post-starburst galaxies, it re- 
mains plausible that they represent a transient phase between an 
intense centrally-concentrated burst of star-formation and a more 
quiescent stage of star-formation activity. 




0.1 
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Fig. 7. The properties of asymmetries and concentrations for 
each class of galaxies (post-starburst with triangles, violet- 
colour coded; quiescent galaxies with squares, orange colour- 
coded; and star-forming galaxies with circles, blue-colour 
coded). In both panels, the vertical dotted lines show the median 
(dashed) and extreme quartiles (dotted) for the lOk-parent sam- 
ple. The extreme quartiles of each distribution are used to show 
the dynamical range of density values in each population. The 
typical error obtained with resampling technique is well within 
the symbols used to plot the median value. 



3.4. Environmental properties 

In this section, we explore the typical environment of post- 
starburst galaxies, considering for example whether they live 
preferentially in overdense regions, and on which physical scales 
(from group- to cluster-scales). We also investigate the relative 
numbers of group galaxies of particular morphologies. This as- 
pect of the post-starburst galaxy population, i.e., the study of 
their environment, is of particular interest insofar as it repre- 
sents an attempt to shed light on the physical mechanisms re- 
sponsible for quenching the star formation in this class of galax- 
ies. Apart from the unprecedented quality and completeness of 
this data set, the zCOSMOS survey provides an ideal sample for 
studying the environmental properties because it is designed to 
detect structures down to the physical scales of galaxy groups 
with a velocity uncert ainty of approxim ately ~ 100 kms -1 , or 

0. 00036 x d + z ) (see lLiilv et al.ll2009h . 

We used the th ree-dimensional density estimators derived by 
iKovac et all d2009h with the ZA DE (Zurich Adaptiv e Density 
Estimator) algorithm. We refer to IKovac etldl ([2009) for a full 
description of the methodology for which we provide only a 
brief outline. Given the current spectroscopic sampling rate 
(~ 30%) of the 10k sample, we attempted to study the entire 
volume targeted by zCOSMOS by adding the 30k photomet- 
ric redshifts of all COSMOS galaxies to the presently available 
10k spectroscopic redshifts. In this way using the reconstructed 
three-dimensional galaxy distribution, we used the density es- 
timates based on a volume-limited catalog after application of 
an adaptive kernel (based in turn on nearest neighbours). This 
choice allowed us to account for the redshift galaxy evolution 
over the redshift range covered in this analysis. We used the 
nearest neighbour to determine the projected local galaxy den- 
sity. The environment of each galaxy relative to the mean density 
at a given redshift (< p(z) >) was expressed as an overdensity 6, 
where 1+5 - p/ < p(z) >. The typical errors associated with this 
estimator are of th e order of log(l +6) = 0.1 -0.15 (see Fig. 4 of 
IKovac et al.l2009t) . As tracers, we used galaxies with magnitudes 
brighter than Mb < (20.5 - z) over the redshifts range sampled 
in this analysis (0.48-1). The typical comoving scales sampled 
range between ~ 2 Mpc (when we use the 5 th nearest neighbour 
estimator) and ~ 8 Mpc (using the 20 th nearest neighbour esti- 
mator). Several estimates of the overdensity of each galaxy are 
available, i.e., weighting by either stellar mass or luminosity. In 
this work, however, we used the number-weighted estimators, 

1. e., no further weighting was applied to our environmental esti- 
mators. 

To investigate the role of the environment on the physical 
scales of groups of gala xies, we used the z COSMOS optically- 
selected group catalog dKnobel et alJ l2009). Both the algorithm 
for group detection and the group catalogues were presented in 
Kno bel et alj ((2009). We describe briefly the main points of the 
algorithm and refer the reader for an exhaustive description to 
the respective paper. Two combined group-finding algorithms 
(the friend-of-friend method and the Voronoi tesselation den- 
sity estimator) were implemented to obtain the final catalogue 
of galaxy groups. This combined approach, tested with Mock 
catalogues, lowers the impact of systematic effects, such as spu- 
rious detections and missing groups, group fragmentation, or 
overmerging. Properties such as velocity dispersion, dynamical 
mass, and group richness, were available for each group with an 
associated error of the order of 20%. The bright magnitude limit 
(at Iab < 22.5) of the survey explains in part the relatively small 
number of detected group members (the majority of them have 
2-4 individual members). The final catalogue lists a total of ap- 
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proximately 800 groups with a completeness ranging between 
70% to 90% (see Fig. 8 of lKnobel et ai1l2009l) and a fraction of 
20% being interlopers. 

There are two peculiarities of the zCOSMOS density field 
(cf. Fig. 18 of lKnobel etail ([2009)) that may affect our group- 
based zCOSMOS results. The first peculiar feature in the group 
reconstruction is a lack of group galaxies in the 10k sample com- 
pared to Mocks at z ~ 0.55. The second peculiarity is an excess 
of group galaxies at z ~ 0.9 that is not reproduced by any of the 
Mock catalogues. These peculiarities might be related to a sig- 
nificant underdensity and overdensity observed at z ~ 0.55 and 
~ 0.9, respe ctively, in the rec onstructed three-dimensional den- 
sity field bv lKovac etal] d2009t) . Finally, we note that at z ~ 0.8, 
the fraction of 10k galaxies in groups is on average about 15%, 
and never exceeds 30%. 

We now study the properties of post-starburst candidate 
galaxies in different environments, and compare their properties 
with those of strong EW[OlI] emitters, quiescent galaxies, and 
the lOk-parent sample. The environmental properties estimated 
using the 5"', 10* and 20 th nearest neighbour for zCOSMOS 
post-starburst galaxies and the comparison samples are shown 
in Fig. [8] The typical physical scales for which our analysis was 
performed, sampled by the three estimators at z e [0.48 - 1.0], 
are ~ 3.5, 5.5, and 8 Mpc, respectively. 

In the left panel of Fig. [8] we plot the distribution of the 
log(l+<5 5 //,), log(l+<J 1 o<i) 1 an d log(l+^20">) (from the left to the 
right panel) for post-starburst galaxies, and for quiescent and 
star-forming galaxies (from bottom to the top panel). The dis- 
tribution of the lOk-parent sample is overplotted for each class 
of galaxy (dashed line). The black circles indicate the respective 
ratio of each class of galaxies to the lOk-parent population as a 
function of the three environmental estimators computed in three 
equally-populated bins of galaxies. 

The bottom panel of Fig. [8] shows a deficit of post-starburst 
galaxies with low values of overdensities (2% with log(l+<5)< 
0.1) that increases to ~ 5% for higher overdensities. Their dis- 
tribution is shifted to higher overdensities than that of the 10k- 
parent population. This trend is supported by the low probabil- 
ity (at a significance level of more than 3cr) that the two dis- 
tributions are drawn from the same lOk-parent population. As 
expected, the fraction of quiescent galaxies also increases at 
higher overdensity values, illustrating their preference to reside 
in denser environments than the global population of galaxies 
(Fig. [8] middle panel). Using the spectroscopic definition de- 
scribed in Sect. 12.21 the fraction of quiescent galaxies doubles 
from 15% in underdense regions (log(l+<S)< 0.1) to 30% in 
overdense regions (> 0.5). The top panel of Fig. [8] shows the 
frequency of highly star-forming galaxies relative to the total 
number of the lOk-parent ga laxies in different de nsity regions. 
As demonstrated in detail by [Cucciati et al.l d2009l) . but see also 
llovino et alj (12009) for the complex interplay between groups 
and galaxy properties, there is no significant correlation between 
the star-forming galaxies in zCOSMOS and their environment- 
It is useful to no t e that the definition of star-forming galaxies 
in ICucciati et al.l ((2009) does not differ significantly from the 
o ne adopted in this wo rk. A star-forming galaxy was defined 
in lCucciati et alJ (12009) to be an object with a SFR higher than 
~ ^Aloyr"" 1 , calibrat ed using the [Oil] luminosity following 
Moustaka s et alj d2006l) . 



2 The possible trend of a higher fraction of star-forming galaxies be- 
ing in underdense environ ment when using the m ass-weighted overden- 
sities has been explored in Cucciati et al. (2009), but it disappears when 
using number densities (as employed here). 



Globally, the distributions of both star-forming and quiescent 
galaxies have a low probability (at > 3cr) of being drawn from 
the same lOk-parent population. 

In the right panel of Fig. [8] we plot the median values of the 
three estimators for all the samples (post-starburst with violet- 
coded triangles, quiescent with orange squares, and star-forming 
galaxies with blue circle symbols). The first (25'' 1 ) and fourth 
(75' /! ) quartiles of each distribution are plotted to show the range 
in density of each population. By dashed, vertical lines, we in- 
dicate the corresponding statistical values of the lOk-parent cat- 
alogue. The typical error bar for environmental quantities in our 
galaxy samples is obtained by a resampling technique and plot- 
ted in the bottom right panel of Fig. [8] The representation of 
median values and the distribution of the central quartiles con- 
firm also that quiescent and post-starburst galaxies reside in an 
overdense environment with respect to lOk-parent galaxies. We 
note that this environmental dependence is observed on different 
physical scales for both galaxy types with a statistically signif- 
icance of more than 3<x level, and that in addition this depen- 
dence in quiescent galaxies is stronger on progressively smaller 
physical scales. On the other hand, star-forming galaxies do not 
appear to reside in significantly different environments from the 
lOk-parent galaxies on all physical scales (from 8 Mpc down to 
~2Mpc). 

We are now able to investigate the impact of the environ- 
ment probed on group physical scales, on the properties of post- 
starburst g alaxies. Using the o ptically selected sample of group 
galaxies bv lKnobel et al.l(l2009l) . we measured and compared the 
fraction of galaxies in each class identified as group members. In 
computing these fractions, we accounted for the redshift depen- 
dence of the group richness by using an effective richness (see 
for details Knobel et al. 2009). 

Using the 10k catalogue of groups with a group defined 
to have either two or more than three members, post-starburst 
galaxies do not appear to exist preferentially in a group en- 
vironment; the fraction of post-starburst galaxies in groups is 
comparable to that of star-forming galaxies (~ 20%) and lower 
than that of quiescent galaxies (using MCOR selection criteria). 
Using the more conservative criterion on mass completeness, 
i.e., selecting only galaxies more massive than log(A1/A1o)> 
10.8, we observe that the fractions of star-forming and post- 
starburst galaxies in groups increase, and become similar to that 
of quiescent galaxies (~ 30%). 

We note that the majority of both post-starburst and quies- 
cent galaxies residing in groups have an early-type morphology 
(~ 71% and ~ 85%, respectively), compared to ~ 17% for the 
star-forming galaxies in sample MCOR. However, for smaller 
numbers of objects, similar values are found for sample MST, 
which support these findings. 

Based on these results, we conclude that although post- 
starburst galaxies are found preferentially in overdense envi- 
ronments on typical scales of clusters and above, they also ex- 
ist in groups and in the field. Given the results for the typical 
k + a galaxy environment using nearest neighbours estimators 
and group catalogues, we found that on average lc + a galax- 
ies prefer physical scales of ~ 2-8 Mpc, i.e., overdense en- 
vironments similar to that of quiescent galaxies. We also iden- 
tified post-starburst galaxies in groups (thus on physical scales 
smaller than Mpc), although they tend to share a group environ- 
ment more similar to that of star-forming galaxies. Given that 
during the short interval of the post-starburst phase (w 1 Gyr), 
a galaxy cannot be transported from an under- to an overdense 
environment, our results imply that the k + a spectral features are 
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Fig. 8. The environmental properties of the three classes of objects using different numbers of neig hbours (5'\ 10"\ and 20 th nearest 
neighbours, or NN). In the left and right panels, we plot the distribution and statistics of the (1+(S 5 »), (I+Siqh,), and (l+£ 2 o») for 
each class of galaxies colour coded according to their classification (post-starburst with violet-colour coded, triangles; quiescent 
galaxies orange colour-coded, squares; and star-forming galaxies blue-colour coded, circles). In both panels, the vertical dotted 
lines (gray colour) show the median and extreme quartiles for the lOk-parent sample. The black circles indicate the ratio of each 
class of galaxies with a particular density to the lOk-parent population computed in three equally-populated bins of galaxies. In the 
right panel, the extreme quartiles of each distribution are used to show the range of density values in each population. The associated 
error bars obtained with resampling technique are plotted in the bottom, right panel. 



produced by a series of mechanisms among the several proposed 
in the past and not necessarily related only to environment. 

Similar results were obtained by Poggianti et alj d2009l) at 
Z = 0.4 - 0.8 using the ESO Distant Cluster Survey, who ro- 
bustly proved the higher incidence of k + a galaxies in high den- 
sity environments (in clusters and in groups with a low fraction 
of [Oil] emitters) than in the field, and pointed out that several 
different processes can produce a k + a spectrum. In contrast, 
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in their study of DEEP2 data found a higher 
fraction of post-starburst galaxies in field rather than group en- 
vironments at lcr significance level. We also note that galax- 
ies in our groups may have a wide range of physical proper- 
ties. This scatter in physical properties reflects t he environmen- 
tal de pendence of post-starburst galaxies (e.g., Poggiant i et alj 
l2006h . As shown in detail by Poggian ti et alj d2009b . k + a galax- 
ies in groups with a low (or high) fraction of [Oil] emitters 
have higher (lower) dependencies on environment. Finally, both 
merging or/and accretion events, and cluster-related mechanisms 
such as ram-pressure gas stripping, harassment, or strangulation, 
are efficient in suppressing star-formation activity in a galaxy. 
Multiple possible origins of k + a galaxies have been proposed 
and we indeed confirm that a wide range of origins is likely. 



3.5. Contribution to the mass assembly and spectral 
properties 

The incidence of post-starburst galaxies in the general galaxy 
population varies between authors, depending on the envi- 
ronment studied, redshifts, selection methods, and properties 



of the various surveys. In cluster environments and at inter- 
mediate redshift, the fraction of post- starburst galaxies varies 
from high values, e . g.. ~ 20 - 25% (iDressler & Gunnl Il983| 
Dres sier et all 119991; IDressler & Gunnl Il992t iTran et all £0031 
12004 120071: IPoggianti et al.ll 19991) . to intermediate values (e.g., 
10%. IPoggianti et al J l2009|). to very low values (~ 0.2% by 
iBalogh et al.ldl999t) ; lYan et all (l2009t) ) both in the field and clus- 
ter environment at intermedia te redshift. At lower redshifts, post- 
starburst galaxies are rare dZabludof f et al.1 119961: iBlake et all 
l2004t iGotol [20071) and reside mostly in the field. It is clearly 
difficult to compare these Ic + a galaxy results for parent samples 
selected from different environments and cosmic times and com- 
piled with different selection criteria. We therefore established 
our own internal means of assessing the contribution of k + a 
galaxies to galaxy stellar mass assembly and galaxy evolution in 
general. 

We found a number density of post-starburst galaxies of 
~ 4.6 + 0.3 x 10~ 5 Mpc 3 and a mass density of ~ 2.6 + 0.2 X 
10 5 M Q Mpc 3 over the redshift range 0.48 < z < 1.2. The er- 
rors were estimated using resampling techniques: we estimate 
the sample statistics f„ using all n galaxies, then remove one 
point at a time to estimate f n -i,i- We found no statistically sig- 
nificant evolution with redshift in either the number or mass den- 
sities of post-starburst galaxies, but given the limited number of 
galaxies in each redshift bin (~ 10) a more solid conclusion will 
be drawn at the completion of the survey. 

To establish the impact of the k + a galaxy population on 
galaxy evolution, we estimated the total mass flux that enters the 
red sequence after a (< 1 Gyr) quenching of the star-formation 
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activity. We followed the method of Mart in et~ai1 (120071) who 
introduced this parameter in their study of SDSS/GALEX red 
and blue galaxies defined with NUV-r colours, which was sub- 
sequently used in other works, e.g., by lArnouts et al.l d2007l) on 
a SWIRE- VVDS-CFHTLS data set. We defined the flow in stel- 
lar mass by post-starburst galaxies that contribute to the red- 
sequence mass assembly as 

1 ^ M* 

Pk+a->RedSeq = " / . 7} ; (1) 

h+a Wj.V max j 

where M* is the stellar mass in solar masses of k + a galax- 
ies in the sample MCOR, w is the selection function (wssr 
and wjsr as introduced in Sect. l2.51 >. and V max is the maxi- 
mum volume in which each galaxy of a given Iab magnitude 
is still observable. We assume that the k + a spectroscopic fea- 
tures are detectable for tu+ a = 0.35 Gyr, or 1 Gyr. Our k + a 
galaxies correspond to a mass flux entering the red-sequence of 
Pk+a^RedSeq = 7.8 x 10~ 4 + 1 .2 M Q Mpc~ 3 yr~' (or a contribution 
of ~ 8% to the growth rate of the red-sequence estim ated to be 
AM/ At ~ 10" 2 Ms Mpc 3 yr _I bv lPozzetti et alj d2009l) and as- 
suming tk+a = 0.35 Gyr, or ~ 3% with %+ a = 1 Gyr). Of course, 
these values are in agreement with the assumption that galax- 
ies remain on the red-sequence and a ny pro cess of re-juvenation 
such as those proposed bv lHasingerl d2008|). are at work. 

Results from the GDDS obtained bv lLe Borgne et alJ J2006) 
showed that galaxies with strong H5 absorption (H6 strong 
galaxies, or HDSs) decrease with cosmic time (from about 50% 
at z ~ 1.2 to a few percent today). However, they detected only 
one or two post-starburst galaxies among the 25 HDSs (depend- 
ing on how their error bars were treated). A strong evolution 
in the number de nsity of post-starburst galaxies was found by 
IWild et all d2009t) (a factor of 20 highe r at z = 0.7 than in the 
local Universe bv lOuintero et al.l (2004)), but on the basis of 5 
post-starburst galaxies selected only with a quite conservative 
definition i n the VVDS-F02 catalogue. Based on that sample, 
IWild et al.l (120091) estimated that post-starburst galaxies can ac- 
count for ~ 40% of the growth in the re d sequence in the field at 
Z = 0.7. Instead. iDe Lucia et al.l (|2009) estimated that the post- 
starburst phase is not the dominant channel that moves galaxies 
from the blue star forming cloud to the red sequence in clusters 
at z = 0.5 -0.8. 

We now investigate the possible evolution in the spectral 
properties of our post-starburst galaxies. We use spectra of post- 
starburst galaxies (in sample MCOR) at high redshift (z e 
[0.75 - 1.2]) and low redshift (z e [0.48 - 0.75]). In the three 
panels of Fig.|9l we show from top to bottom the composite 
spectra of post-starburst galaxies (in sample MCOR) at high red- 
shift, low redshift, and the difference between them. The redshift 
bins are chosen to include a comparable number of post-starburst 
galaxies (40 spectra at low z, and 28 at high z), although the spec- 
tra are weighted using our selection function before the combi- 
nation. The median redshifts of the two composite spectra are 
Zmed = 0.659 and Zmed = 0.895 at low and high redshifts, respec- 
tively. We shifted our spectra to the galaxy rest-frame defined by 
their spectroscopic redshifts, and normalized them to the com- 
mon region around 4100^-700/1 sampled in each spectrum. To 
each spectrum, we applied the weighting approach described in 
Sect. [23] obtaining a composite spectrum with a median combi- 
nation. 

We do not observe any significant difference between the 
spectral properties of post-starburst galaxies at Zmed = 0.659 and 
Zmed = 0.895. The continuum of the difference spectrum is cen- 
tred on zero at all wavelengths (apart from at very blue A, where 



the difference reaches 5%). It is, however, true that there is not 
a very large difference in redshift between these two spectra, so 
negligible evolution given the uncertainty levels, may not be a 
surprising result. 
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Fig. 9. Composite spectra of post-starburst galaxies at z e [0.75 - 
1.2] and z e [0.48 - 0.75] (and the non-statistically difference 
between them in the bottom panel). 



4. Discussion 

A number of studies have explored the physical properties and 
the relative importance of k+a galaxies, i.e., galaxies with strong 
Balmer absorption lines (which are indicative of an intense star- 
formation epoch in the past billion years) and an absence of 
emission lines (a signature of lack of ongoing star formation). 
Even considering the different numbers cited in the literature 
about the incidence of k + a galaxies among the entire galaxy 
population, and their apparent or true evolution in number den- 
sity (which has not been confirmed in this work), they represent a 
minor constituent of the galaxy population in the Universe in the 
past 8 Gyr. These systems are physically interpreted as galaxies 
that have had their star formation abruptly quenched in the past 
billion years and subsequently migrated onto the red sequence. 
Thus, the importance of this galaxy population depends on its 
role in creating the red sequence, or in stellar mass assembly in 
general. Other important issues are the mechanism/s responsible 
for the abrupt quenching of star formation in their progenitors 
and the true nature of their descendants. 

We have completed an analysis of high quality spectral data 
provided by zCOSMOS and ancillary data available for the 
COSMOS field of the first mass-selected sample of 69 k + a 
galaxies in a wide range of environments at intermediate redshift 
(0.48 <z< 1.2). 

Alternative explanations of why k + a galaxies appear to be 
transient objects have also been considered. Using 1 .4 GHz ra- 
dio and X-ray data, we have confirmed that our k + a galax- 
ies are unlikely to be strongly star-forming galaxies affected 
significantly by dust extinction and only a fraction of them 
may be obscured AGNs. We have also considered another con- 
cern that Balmer line strengths might be generated by some 
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other hot population tha n the main-sequen ce turn-off stars, 
e.g., blue stragglers (e.g.. lRosdll985L Il994t) or blue horizon- 
tal branch stars (e.g..lBurstein et alj|1984t iMaraston & Thomasl 
l2000HTrager et al.ll2005h . 

In this work, we adopted a conventional definition in se- 
lecting k + a galaxies. Although Ha is a robust indicator of 
star formation, it falls out of the optical window at interme- 
diate redshifts and bluer diagnostics are used in assessing the 
star-formation a ctivity of the post-starburst galaxy population. 
lYan et al.l d2006l) reported on the incompleteness introduced in 
using [Oil] to identify k + a galaxies, a criterion that could un- 
fairly eliminate post-starburst galaxies with AGN activity. This 
is because the [Oil] emission line is sensitive to both the star 
formation and AGN activity in a galaxy, and by including galax- 
ies on the basis of the lack of an [Oil] emission line, we are a 
priori also removing those with an AGN contribution. However, 
we followed the classical approac h of using t he [O il] diagnostic 
instead of the H/3 line proposed by Yan et al.l d2009l) . because H/3 
is a far weaker line than [Oil] which can increase the uncertainty 
in e.g., the stellar continuum fit, and is affected by both emission 
and absorption. Thus, in this paper we adopt the conventional 
definition that involves both the [Oil] emission line and H6a, 
to maximize the comparison with the majority of other works. 
Any impact of the k + a galaxy population in a general con- 
text of galaxy evolution may be slightly greater than quantified 
here. Our incompleteness should however be minimal because 
the duty cycle (~ 10 8 yr, e.g., Condon et al. 1992) of a possible 
AGN in our k+a galaxies represents a short period in the lifetime 
of a galaxy. 

We note that any comparison with external results is possi- 
ble only qualitatively, while a more objective analysis is feasible 
with internal comparisons within a given catalogue as pointed 
out by several authors. We followed this approach in our work 
by objectively comparing our results on k + a galaxies with star- 
forming and quiescent galaxy mass-selected samples defined 
using spectroscopic criteria and with lOk-parent mass-selected 
galaxies, all extracted from the zCOSMOS bright program. The 
results obtained in this work are robust and corrected for effects 
of incompleteness. 

We explored the physical properties of k+a galaxies in terms 
of stellar mass, mass-to-light ratio, colour, morphology, and en- 
vironment. We found that our k + a galaxies are as massive as 
quiescent galaxies, but have relatively lower mass-to-light ratio 
and bluer colours, which ensure that they are potential contrib- 
utors to the red sequence and ideal occupants of the green val- 
ley, i.e., the underpopulated region in-between the distribution of 
star-forming/blue galaxies and quiescent/red galaxies. We con- 
firmed that morphologically zCOSMOS k + a galaxies are a het- 
erogeneous population of bulge- and disc-dominated galaxies, 
which become predominantly more bulge-dominated at increas- 
ing stellar mass. Despite the high incidence of asymmetries in- 
dicative of recent interaction events, our morphological results 
have demonstrated that not all progenitors of our k + a galaxies 
are spheroidals as expected if a major merger had destroyed the 
disc. 

Our finding that 50% of our k + a population are passive spi- 
ral galaxies (but depending on stellar mass) and the knowledge 
that in distant clusters, k + a galaxies ar e found to have pre- 
dominantly spiral m orphologies (e.g., by iDressler et al.l ll999: 
Poggi anti et al.lll999t) is consistent with the timescale for mor- 
phological transformation being longer than that for the trans- 
formation in spectral properties. However, because of their high 
stellar masses and the large spread in their colours in the green 
valley, star formation in at least a fraction of k + a galaxies may 



be periodically triggered by the rapid accretion of outer fresh 
material ( as e.g ., proposed by the re-juvenation hypothesis of 
lHasin ger (2008)). Given that k + a galaxies are found in a wide 
range of environments (from field to cluster scales), the mecha- 
nisms quenching star formation in k + a galaxy progenitors and 
their descendants might both be diverse. Summarizing, the re- 
sults obtained in the present work provides good support to the 
common interpretation for k + a galaxies that they may represent 
a transient phase between star-forming and a quiescent stage of 
star-formation activity, but their subsequent evolution may be to 
a wide range of galaxy classes by means of many different evo- 
lutionary paths. 

It is important to determine the way in which star formation 
is controlled and truncated in a galaxy. Star formation can ex- 
perience a gradual decrease as would be consistent for example 
with the smooth decline in redshift (z < 1) of the star forma- 
tion rate, or star formation can be suddenly and violently ex- 
tinguished. Given that the spectra of k + a galaxies imply that 
they have had their star formation truncated suddenly in the re- 
cent past, measuring their number and mass density can help us 
to quantify the significance of such phenomenon in suppress- 
ing the global star formation rate. K + a galaxies have abruptly 
quenched their star formation activity on a timescale shorter than 
a billion years, therefore their contribution represents a lower 
limit to the prevalence of such sudden quenching. The decline in 
the star formation activity in galaxies is understood to have oc- 
curred gradually with cosmic time in the past 8 Gyr, although its 
depe ndence on the stella r m ass remains controversial , e.g., com- 
pare |Zheng!eLaD (120071) and iDrorv & Alvar ez (2008). We com- 
puted the contribution of k + a galaxies to the mass assembly 
of the red-sequence in the conservative hypothesis that they re- 
main on the red-sequence. This contribution to the growth rate of 
the red-sequence has been found to be in the range between 3% 
and 8% depending on the assumed timescale of the characteristic 
detectability of k + a spectroscopic features. Considering a num- 
ber of possible underestimation effects such as e.g., the removal 
of k + a galaxies with obscured AGNs and the non-detection 
of k + a objects because of the insufficient signal-to-noise ra- 
tio of the spectra, we conclude that although the k + a galaxy 
population is not the primary contributor to the red sequence, it 
might not be negligible in the general context of galaxy evolu- 
tion. Although sudden truncations of the star formation activity 
appear not to be more prevalent than smoother declines in the 
redshift interval 0.48 < z < 1.2, they are probably still impor- 
tant. In addition, a lack of sufficient temporal resolution in the 
star formation history cannot prevent us from speculating that 
the star formation might occur in episodes of moderate intensity 
separated by short quiescent phases as observed in nearby dwarf 
galaxies in whi ch the resolved stell ar populations can be studied 
in great detail dTolstov et al.l l2009h . In this scenario, the k + a 
spectral features are common in the life of each galaxy, but ob- 
served only in a small fraction of the entire galaxy population 
given their short lifetime. 



5. Conclusions 

We have explored the physical properties of the first mass- 
selected sample of k + a galaxies at intermediate redshift (z = 
0.48 - 1.2) in a wide range of environments. We selected this 
galaxy population, interpreted to represent a link between star- 
forming/blue galaxies and quiescent/red galaxies, from the spec- 
troscopic program zCOSMOS related to the largest multiwave- 
length survey completed to date, COSMOS. 
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Our results can be summarized as follows: 



1 . K + a galaxies occupy the brightest tail of the luminosity dis- 
tribution. They can be as massive as typical quiescent galax- 
ies and populate the so-called "green valley" in the colour 
versus luminosity (or stellar mass) distribution. 

2. A fraction (< 8%) of these galaxies have radio and/or X- 
ray counterparts implying an upper limit to the SFR of ~ 
8A4 yr~' . These results suggest a possible contribution from 
obscured AGNs. 

3. K + a galaxies morphologically represent a heterogeneous 
population with a similar incidence of bulge-dominated and 
disky galaxies, independent of the strength of their HS ab- 
sorption line, but dependent on stellar mass in a way rem- 
iniscent of the well-known mass-morphology relation (i.e., 
the most massive post-starburst galaxies are predominantly 
bulge-dominated systems). We reported a high incidence of 
asymmetries in HST/ACS images, supporting the hypothe- 
sis that these galaxies have experienced interaction/merger 
processes in their recent past, although we are unable to as- 
certain whether these processes are entirely responsible for 
the quenching of the galaxy star-formation activity. 

4. K+a galaxies reside preferentially in quite rich environments 
on physical scales of ~ 2 - 8 Mpc and have morphological 
early-to-late type ratios in groups that are similar to those 
of quiescent galaxies, although we also found k + a galax- 
ies in underdense regions. We concluded that several differ- 
ent mechanisms, not necessarily related to a specific envi- 
ronment, quench star-formation activity in galaxies on short 
timescales (< 1 Gyr). 

5. We measured a number density of k + a galaxies of ~ 
4.6 ± 0.3 x 1(T 5 Mpc~ 3 and a mass density of ~ 2.6 ± 0.2 x 
10 5 Mo Mpc 3 over the redshift range 0.48 < z < 1.2. With 
the current data set, we did not find statistical evolution in 
either any of the aforementioned quantities in the redshift in- 
terval explored, nor b y comparison with the local values of 
lOuintero et all £2004).' We were also unable to determine any 
evolution in the spectral properties between 0.48 < z < 0.75 
and 0.75 < z < 1.2. 

6. We studied the contribution of k + a galaxies to the mass as- 
sembly of the red-sequence. This contribution can be as high 
as ~ 8% {pk^RedSeq = 7.8 x 10- 4 ± 1.2 M & Mpc" 3 yr 1 ) 
given a growth rate in the red-sequence of dM/dt ~ 
10~ 2 M e Mpc~ 3 yr~' (as estimated by Pozzetti et al. 2009 
for the same data) and assuming a detectability of k + a spec- 
troscopic features of tk+ a - 0.35 Gyr. Although k + a galax- 
ies may not be the primary source of growth to the red se- 
quence, they are likely to be important to galaxy evolution. 
Furthermore, sudden truncations in star formation do not ap- 
pear to be the primary means of quenching the star formation 
activity compared to more gradual means, but appear never- 
theless to be important to the decline in the global star forma- 
tion activity. On the other hand, lacking sufficient temporal 
resolution, we cannot exclude the scenario where the sudden 
quenching of star formation and the regime in which the star 
formation occurs in episodes of moderate intensity separated 
by short quiescent phases are the rule in general galaxy evo- 
lution. It is worth noting indeed that in nearby dwarf galaxies 
when the resolved stellar populations are studied in great de- 
tail, the star formation history of these galaxi es appears to 
consi st of repeated periods of star formation (Tolsto y" ;t al.l 
2009). More solid results on this and other subjects will be 
obtained using the larger 20k zCOSMOS sample. 



6. Acknowledgments 

This work has been supported in part by the grant 
ASI/COFIS/WP3110 1/026/07/0. We thank the referee for pro- 
viding constructive comments. 

References 

Arnouts, S., Walcher, C. J., Le Fevre, O., et al. 2007, A&A, 476, 137 
Balogh, M. L., Miller, C, Nichol, R., Zabludoff, A., & Goto, T. 2005, MNRAS, 
360, 587 

Balogh, M. L. & Morris, S. L. 2000, MNRAS, 318, 703 

Balogh, M. L., Morris, S. L., Yee, H. K. C, Carlberg, R. G., & Ellingson, E. 
1999, ApJ, 527, 54 

Barger, A. J., Aragon-Salamanca, A., Ellis, R. S., et al. 1996, MNRAS, 279, 1 
Barnes, J. E. & Hernquist, L. 1992, ARA&A, 30, 705 

Bell, E. E, Mcintosh, D. H., Katz, N, & Weinberg, M. D. 2003, ApJS, 149, 289 
Blake, C, Pracy, M. B., Couch, W. J., et al. 2004, MNRAS, 355, 713 
Bolzonella, M., Kovac, K., Pozzetti, L., et al. 2009, ArXiv e-prints 0907.0013 
Bondi, M., Ciliegi, P., Schinnerer, E., et al. 2008, ApJ, 681, 1129 
Borch, A., Meisenheimer, K., Bell, E. E, et al. 2006, A&A, 453, 869 
Bower, R. G., Kodama, T., & Terlevich, A. 1998, MNRAS, 299, 1193 
Brinchmann, J., Chariot, S., White, S. D. M., et al. 2004, MNRAS, 351, 1151 
Brinchmann, J. & Ellis, R. S. 2000, ApJ, 536, L77 
Brusa et al. 2009, in prep. 

Bundy, K., Ellis, R. S., Conselice, C. J., et al. 2006, ApJ, 651, 120 

Burstein, D., Faber, S. M., Gaskell, C. M., & Krumm, N. 1984, ApJ, 287, 586 

Caldwell, N. & Rose, J. A. 1997, AJ, 113, 492 

Caldwell, N., Rose, J. A., & Dendy, K. 1999, AJ, 117, 140 

Cappelluti et al. 2009, in prep. 

Cayatte, V., van Gorkom, J. H., Balkowski, C, & Kotanyi, C. 1990, AJ, 100, 604 
Chabrier, G. 2003, PASP, 1 15, 763 

Chang, T.-C, van Gorkom, J. H., Zabludoff, A. I., Zaritsky, D., & Mihos, J. C. 

2001, AJ, 121, 1965 
Chariot, S. & Fall, S. M. 2000, ApJ, 539, 718 
Condon, J. J. 1992, ARA&A, 30, 575 

Conselice, C. J., Bershady, M. A., & Jangren, A. 2000, ApJ, 529, 886 

Couch, W. J., Ellis, R. S., Sharpies, R. M., & Smail, I. 1994, ApJ, 430, 121 

Couch, W. J. & Sharpies, R. M. 1987, MNRAS, 229, 423 

Cucciati et al. 2009, submitted to A&A 

De Lucia, G. & Blaizot, J. 2007, MNRAS, 375, 2 

De Lucia, G., Poggianti, B. M., Halliday, C, et al. 2009, MNRAS, 1379 
Dewangan, G. C, Singh, K. P., Mayya, Y. D., & Anupama, G. C. 2000, MNRAS, 
318, 309 

Dressier, A. & Gunn, J. E. 1983, ApJ, 270, 7 
Dressier, A. & Gunn, J. E. 1992, ApJS, 78, 1 

Dressier, A., Oemler, A. J., Sparks, W. B., & Lucas, R. A. 1994, ApJ, 435, L23 

Dressier, A., Smail, I., Poggianti, B. M., et al. 1999, ApJS, 122, 51 

Drory, N. & Alvarez, M. 2008, ApJ, 680, 41 

Elvis, M., Civano, E, Vignali, C, et al. 2009, ApJS, 184, 158 

Fontana, A., Pozzetti, L., Donnarumma, I., et al. 2004, A&A, 424, 23 

Franx, M. 1993, ApJ, 407, L5 

Franzetti, P., Scodeggio, M., Garilli, B., et al. 2007, A&A, 465, 71 1 
Georgakakis, A., Nandra, K., Yan, R., et al. 2008, MNRAS, 385, 2049 
Giovanelli, R. & Haynes, M. P. 1985, ApJ, 292, 404 
Goto, T. 2005, MNRAS, 357, 937 
Goto, T. 2007, MNRAS, 377, 1222 

Goto, T. & SDSS collaboration. 2003, in Bulletin of the American Astronomical 

Society, Vol. 35, Bulletin of the American Astronomical Society, 1249 
Gunn, J. E. & Gott, J. R. I. 1972, ApJ, 176, 1 
Hammer, E, Flores, H., Lilly, S. J., et al. 1997, ApJ, 481, 49 
Hasinger, G. 2008, A&A, 490, 905 

Hopkins, P. F, Bundy, K., Hernquist, L., & Ellis, R. S. 2007, ApJ, 659, 976 
Iovino et al. 2009, ArXiv e-prints 0909.1951 

Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003a, MNRAS, 341, 
33 

Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003b, MNRAS, 341, 
54 

Kaviraj, S., Kirkby, L. A., Silk, J., & Sarzi, M. 2007, MNRAS, 382, 960 
Kim et al. 2009, in prep. 

Knobel, C, Lilly, S. J., Iovino, A., et al. 2009, ApJ, 697, 1842 
Koekemoer, A. M., Aussel, H, Calzetti, D., et al. 2007, ApJS, 172, 196 
Koopmann, R. A. & Kenney, J. D. P. 2004, ApJ, 613, 866 
Kovac, K, Lilly, S. J., Cucciati, O., et al. 2009, ArXiv e-prints 0903.3409 
Lamareille, E, Contini, T., Le Borgne, J.-E, et al. 2006, A&A, 448, 893 
Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 
Le Borgne, D., Abraham, R., Daniel, K., et al. 2006, ApJ, 642, 48 



D. Vergani et al.: Post-starburst galaxies in COSMOS field 



17 



Lilly, S. J., Le Fevre, O., Renzini, A., et al. 2007, ApJS, 172, 70 
Lilly, S. J., LeBrun, V., Maier, C, et al. 2009, ApJS, 184, 218 
Liu, C. T. & Kennicutt, Jr., R. C. 1995, ApJS, 100, 325 
Maraston, C. & Thomas, D. 2000, ApJ, 541, 126 

Martin, D. C, Wyder, T. K., Schiminovich, D., et al. 2007, ApJS, 173, 342 
Meneux, B., Guzzo, L., de la Torre, S., et al. 2009, ArXiv e-prints 0906.1807 
Mignoli, M., Zamorani, G., Scodeggio, M., et al. 2009, A&A, 493, 39 
Miller, N. A. & Owen, F. N. 2002, AJ, 124, 2453 

Miyaji, T., Griffiths, R. E., & C-COSMOS Team. 2008, in AAS/High Energy 

Astrophysics Division, Vol. 10, AAS/High Energy Astrophysics Division 
Moustakas, J., Kennicutt, Jr., R. C, & Tremonti, C. A. 2006, ApJ, 642, 775 
Norton, S. A., Gebhardt, K., Zabludoff, A. I., & Zaritsky, D. 2001, ApJ, 557, 150 
Poggianti, B. M., Aragon-Salamanca, A., Zaritsky, D., et al. 2009, ApJ, 693, 1 12 
Poggianti, B. M., Bressan, A., & Franceschini, A. 2001, ApJ, 550, 195 
Poggianti, B. M., Smail, I., Dressier, A., et al. 1999, ApJ, 518, 576 
Poggianti, B. M., von der Linden, A., De Lucia, G., et al. 2006, ApJ, 642, 188 
Pozzetti, L., Bolzonella, M., Lamareille, F, et al. 2007, A&A, 474, 443 
Pozzetti, L., Bolzonella, M., Zucca, E., et al. 2009, ArXiv e-prints 0907.5416 
Press, W. H., Teukolsky, S. A., Vetterling, W. T, & Flannery, B. P. 
1992, Numerical recipes in FORTRAN. The art of scientific computing 
(Cambridge: University Press, 1992, 2nd ed.) 
Quintero, A. D., Hogg, D. W., Blanton, M. R., et al. 2004, ApJ, 602, 190 
Rose, J. A. 1984, AJ, 89, 1238 
Rose, J. A. 1985, AJ, 90, 1927 
Rose, J. A. 1994, AJ, 107, 206 

Salim, S., Chariot, S., Rich, R. M., et al. 2005, ApJ, 619, L39 
Scarlata, C, Carollo, C. M., Lilly, S., et al. 2007a, ApJS, 172, 406 
Scarlata, C, Carollo, C. M., Lilly, S. J., et al. 2007b, ApJS, 172, 494 
Schinnerer, E., Smolcic, V., Carilli, C. L., et al. 2007, ApJS, 172, 46 
Schmidt, M. 1968, ApJ, 151, 393 

Scodeggio, M., Franzetti, P., Garilli, B., et al. 2005, PASP, 117, 1284 
Scoville, N, Aussel, H., Brusa, M., et al. 2007, ApJS, 172, 1 
Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 776 
Tasca, L. A. M., Kneib, J. -P., Iovino, A., et al. 2009, A&A, 503, 379 
Tolstoy, E., Hill, V., & Tosi, M. 2009, ArXiv e-prints 0904.4505 
Toomre, A. & Toomre, J. 1972, ApJ, 178, 623 

Trager, S. C, Worthey, G., Faber, S. M., & Dressier, A. 2005, MNRAS, 362, 2 
Tran, K.-V. H., Franx, M., Illingworth, G., Kelson, D. D., & van Dokkum, P. 

2003, ApJ, 599, 865 
Tran, K.-V. H., Franx, M., Illingworth, G. D., et al. 2007, ApJ, 661, 750 
Tran, K.-V. H., Franx, M., Illingworth, G. D., et al. 2004, ApJ, 609, 683 
Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898 
Vergani, D., Scodeggio, M., Pozzetti, L., et al. 2008, A&A, 487, 89 
Wild, V, Walcher, C. J., Johansson, P. H., et al. 2009, MNRAS, 395, 144 
Worthey, G. & Ottaviani, D. L. 1997, ApJS, 111, 377 
Yan, R., Newman, J. A., Faber, S. M., et al. 2009, MNRAS, 398, 735 
Yan, R., Newman, J. A., Faber, S. M., et al. 2006, ApJ, 648, 281 
Yang, Y, Zabludoff, A. I., Zaritsky, D., Lauer, T. R., & Mihos, J. C. 2004, ApJ, 

607, 258 

Yang, Y, Zabludoff, A. I., Zaritsky, D., & Mihos, J. C. 2008, ApJ, 688, 945 

Zabludoff, A. I., Zaritsky, D., Lin, H., et al. 1996, ApJ, 466, 104 

Zheng, X. Z., Bell, E. F, Papovich, C, et al. 2007, ApJ, 661, L41 

Zucca, E., Bardelli, S., Bolzonella, M., et al. 2009, ArXiv e-prints 0909.4674 



13 INAF - Osservatorio Astronomico di Torino, 10025 Pino Torinese, 
Italy 

14 Space Telescope Science Institute, 3700 San Martin Drive, 
Baltimore, MD 21218 

15 LBNL & BCCP, University of California, Berkeley, CA 94720, 
USA 

16 Centre de Physique Theorique, Marseille, Marseille, France 

17 Institut d'Astrophysique de Paris, UMR 7095 CNRS, Universite 
Pierre et Marie Curie, 98 bis Boulevard Arago, F-75014 Paris, 
France. 

18 Universitats-Sternwarte, Scheinerstrasse 1, D-81679 Muenchen, 
Germany 

19 Argelander-Institut fur Astronomie, Auf dem Hiigel 71, D-53121 
Bonn, Germany 

20 INAF, Osservatorio di Roma, Monteporzio Catone (RM), Italy 

21 California Institute of Technology, MC 105-24, 1200 East 
California Boulevard, Pasadena, CA 91 125, USA 

22 Institute for Astronomy, University of Hawaii, 2680 Woodlawn 
Drive, Honolulu, HI, 96822 

23 Research Center for Space and Cosmic Evolution, Ehime 
University, Bunkyo-cho, Matsuyama 790-8577, Japan 



1 INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1, I- 
40127, Bologna, e-mail: daniela.vergani@oabo.inaf.it 

2 Universita di Bologna, Dipartimento di Astronomia, Via Ranzani 
1, 1-40127, Bologna, Italy 

3 Institute of Astronomy, Swiss Federal Institute of Technology 
(ETH Honggerberg), CH-8093, Zurich, Switzerland. 

4 Laboratoire d'Astrophysique de Toulouse-Tarbes, Universite de 
Toulouse, CNRS, 14 avenue Edouard Belin, F-31400 Toulouse, 
France 

5 INAF - Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 
1-50125 Firenze, Italy 

6 INAF - IASF Milano, Milan, Italy 

7 INAF - Istituto di Radioastronomia, Via Gobetti 101, 1-40129 
Bologna, Italy 

8 INAF Osservatorio Astronomico di Brera, Milan, Italy 

9 Laboratoire d'Astrophysique de Marseille, Marseille, France 

10 European Southern Observatory, Karl-Schwarzschild-Strasse 2, 
Garching, D-85748, Germany 

1 1 Dipartimento di Astronomia, Universita di Padova, Padova, Italy 

12 Max-Planck- Institut fur extraterrestrische Physik, D-84571 
Garching, Germany 



